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ABSTRACT 

Many massive galaxies at the centers of relaxed galaxy clusters and groups have vast reser¬ 
voirs of cool (^10,000 K) and cold (^100 K) gas. In many low redshift brightest group 
and cluster galaxies this gas is lifted into the hot ISM in filamentary structures, which are 
long lived and are typically not forming stars. Two important questions are how far do these 
reservoirs cool and if cold gas is abundant what is the cause of the low star formation effi¬ 
ciency? Heating and excitation of the filaments from collisions and mixing of hot particles in 
the surrounding X-ray gas describes well the optical and near infra-red line ratios observed 
in the filaments. In this paper we examine the theoretical properties of dense, cold clouds 
emitting in the far infra-red and sub-millimeter through the bright lines of [C ll]A157/im, 
[O I] A63/im and CO, exposed to these energetic ionising particles. While some emission lines 
may be optically thick we find this is not sufficient to model the emission line ratios. Models 
where the filaments are supported by thermal pressure support alone also cannot account for 
the cold gas line ratios but a very modest additional pressure support, either from turbulence 
or magnetic fields can fit the observed [O l]/[C ll] line ratios by decreasing the density of 
the gas. This may also help stabilise the filaments against collapse leading to the low rates of 
star formation. Finally we make predictions for the line ratios expected from cold gas under 
these conditions and present diagnostic diagrams for comparison with further observations. 
We provide our code as an Appendix. 
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1 INTRODUCTION 

The Herschel telescope opened up a wealth of lines in the far 
infra-red (FIR) such as the strong fine structure cooling lines 
of [O l]A63/imand [C ll]A157/ini. These emission lines, which 
are often attributed to the presence of photo-dissociation regions 
(PDRs), excited by far Ultra-Violet (FUV) emission from young 
stars (e.g. |Hollenbach & Tielens|1999| l, have been found in a num¬ 
ber of brightest cluster galaxies (BCGs, e.g. [Edge et aL]|2010a|b[ 
[Mittal et al.|201 1|[20 12 [[Werner et al.|2013|) and in X -ray and radio 
bright giant elliptical (gE) galaxies ([Wemer et al.[[2014| [Guillard[ 

[et al.|2014[ l. In these dense. X-ray bright systems the cooling time 
of the hot gas is short; an important consequence of which is that a 
heat source is required to prevent catastrophic cooling of the intra- 
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cluster medium (ICM) which would lead to unprecedented growth 
of the galaxy. This heating is likely provided by active galactic nu¬ 
cleus (AGN) feedback (for a review see [Fabian[2012[ ). 


Though sample sizes are small, where extended ionised 
emission lines have been observed, Herschel observations of 
[C ll]A157/imhave shown the cold (T<100 K) gas is also ex¬ 
tended and coincident with filaments of ionised gas ( [Werner[ 
[et al.[[20T3l [2014[ ). [Werner et al.[ ( [2014[ ) also show that the 
Ha/[C ll]A157/imratio is similar implying that the ionised and 
cold gas share a common excitation mechanism. Though evidence 
for some star formation is observed in BCGs (e.g. [Hubble & Hu 
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Figure 1. Left: The temperature of the cloud. For the optically-thin case the temperature on this plot is the temperature of e ach ‘cloudlet’ and in the optically- 
thick case it will be the surface temperature of the cloud. The solid black line represents a total pressure of 10®-^ cn^^K (the gas pressure, derived from 
X-ray measurements, of the core of the Perseus cluster) and the hatched area represents pressures from 10^-10^ cnEElK . Middle: Contours of the total gas 
pressure in the simple optically-thin, non-turbulent case with no additional pressure from magnetic fields. Right: Ionisation fraction, Xe, (electron density, rie, 
divided by particle density, n) of the gas for the same model. 
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20121, many of these extended filaments have no obvious star for¬ 
mation despite the large quantities of cool and cold gas. PDRs are 
therefore an unlikely heating mechanism. An important question is 
why are these cold gas rich systems not forming stars? 

The ionised filaments are thought to be long-lived but their 
energy source and the reason for their quiescence has been a puz¬ 
zle. H 2 , CO, dust and polycyclic aromatic hydrocarbon (PAH) fea¬ 
tures have been shown to exist coincident with some ionised fil- 
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Bremer & Baker[[2005[ [Crawforc 
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2011[ [Canning et al.||2013] l. The new information garnered from 
the low-ionisation and molecular emission lines seen by Herschel 
and the high spatial resolution made available by ALMA will be in¬ 
strumental in probing density and temperature sensitive lines with 
which to characterise the properties of the cold gas, and in allow¬ 
ing the morphology and kinematics of this gas to be mapped on the 
same spatial scales as the optically emitting ionised filaments (e.g. 
[McNamara et al.|2()T4l [Russell et al.|2()T4) . 

[Ferland et al.[ ( [2008[ ) and [Ferland et al.[ ( [2009[ ) (hereafter 
F2009) have suggested collisions with the surrounding energetic 
particles as a mechanism for heating the multi-phase filaments. The 
authors show that their model, which allows the filaments to be de¬ 
scribed by cloudlets of varying densities, can successfully repro¬ 
duce the ratios of the near infra-red (NIR) strong H 2 , and optical 
atomic and low ionisation emission lines, which has been a struggle 
for photoionisation models. The model also reproduces the charac¬ 
teristic signatures of relatively high He I and [N l] emission and of 
relatively low [O III] in the extended filaments ([Johnstone & Fabian| 


1988[|Wit & Donahue[[1997| Sabra, Shields & Filippenko[[2000[ 

Hatch et al.[[2005[ Canning et al.[[201l| l, while allowing anoma- 

lously high [Ne III] emission, through charge-transfer reactions. 
However, the models are unable to match recent Herschel obser¬ 
vations of the ratios of [O l]A63/im and [C ll]A157/im in the fila¬ 
ments. It should be noted that the authors stress their model is ap¬ 


plicable only to optically-thin clouds and additionally that the data 
used in F2009 cannot constrain a large reservoir of cold clouds in 
the filaments as the power-law controlling the composition of the 
cloudlets is hinged on emission lines from gas at 10^ — 1000 K. 


Collisional excitations and ionisations of atomic and molecu¬ 
lar gas by high energy particles is probably important, not only in 
the filaments of BCGs, but in a wide range of astrophysical environ¬ 
ments such as the opaque molecular cores of PDRs where heating 
is thought to be dominated by excitations produced by cosmic rays 
(e.g. |Hollenbach & Tielens 1997 1 , in cool and cold gas near AGNs 
(e.g. [Shull & van Steenberg 1985 1 and also in the cool gas filaments 
expelled in the death throes of supernovae where non-thermal elec¬ 
trons are generated by the hard synchrotron photoionising source 
(e.g. [Richardson et al.|2013[ l. 


We have extended the models of F2009 to explore the physical 
and chemical state of the cold (T<100 K) gas with densities rang¬ 
ing from n = 10° — 10° cm“° , and in particular we consider the 
effects of finite column densities of individual cloudlets. Our pri¬ 
mary aim is to explore the discrepancy between the observed and 
predicted line ratios of the FIR fine structure lines of [O l]A63/im 
and [C ll]A157/im in the extended filaments of BCGs. However, 
our models are generally applicable to ionised, neutral and molec¬ 
ular gas exposed to energetic ionising particles. The spectrum from 
collisional heating in cool and cold gas clouds is complicated as it 
depends sensitively on changes in the heating and cooling of the 
gas which in turn depends on much of the physics of the cloud. 


Section reviews briefly the details of the particle heating 
model grid and discusses the physical and chemical properties of 
the grid. Section [4J^ explores the effect of a large quantity of cold 
clouds on the composite model spectrum for BCGs. Section [4~2| dis- 
cusses the effects of optical depth on the predictions for the BCG 
line ratios and Sections |43] and [T^ explore the effect of turbulence 
and magnetic fields. Some words of caution about the models is 
given in Section [43] A discussion and predictions for future obser¬ 
vations is presented in Section]^ 
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2 COLD CLOUDS IN COOL CORES 

As mentioned in Section massive central galaxies in relaxed 
galaxy clusters and groups often host large amounts of cool and 
cold gas (e.g. [Crawford et al.|19^ . The dichotomy between the 
lack of gas observed in central galaxies in non-cool core clus¬ 
ters and its abundance in cool core clusters indicates that cooling 
from the X-ray gas must be important in these systems (e.g. jCav-] 
jagnolo et al.||2008| ). However, why these extended gas reservoirs 
continue to shine has been an issue of contention; optical emis¬ 
sion line ratios are typically ‘LINER-like’ (low ionisation nuclear 
emission-line region). However, it has been shown that many phys¬ 
ical processes can produce similar liner-like ratios (weak shocks 


e.g. Sparks, Macchetto & Golombek||1989| [Farage, McGregor &[ 

[Dopita|2012| hot stars e.g. 

Terlevich & Melnick|1985| conduction 

e.g.[Voit & Donahue[1997| 

Sparks et al.|2012[ mixing with hot gas 


which leads to both collisonal excitation and thermal excitation of 


the gas e.g. F2009; [Fabian et al.|2011| and dissipation from mag¬ 
netic reconnection in magnetically supported filaments e.g. [Chur^ 
[zov, Ruszkowski & Schekochihin[2013] l. 

The best, nearby, example of an extended ionised and molec¬ 
ular web surrounding a BCG is NGC 1275. Some regions of this 
extended web are forming stars but the majority are not ( [Conselice^ 
[Gallagher & Wyse|[2001[ [Canning et ^[2010[ [2014[ l. Obviously, 
where cool gas is close to intense photoionisation from stars, such 
as in the inner regions of some BCGs or in star-forming filaments, 
this photoionisation will contribute significantly to the excitation of 
the gas. We must therefore look at ‘clean’ regions, devoid of obvi¬ 
ous star formation in order to examine the excitation mechanism 
for the filaments prior to their degeneration into stars, or evapora¬ 
tion into the hot ISM/ICM. HST SBC data shows that the SFR in a 
‘typical’ non star forming region, of NGC1275’s filaments, assum¬ 
ing all the UV fiux is due to stars, is <0.001 M© yr 
[et al.|201^ . 

The two questions of basic importance are how much cold gas 
exists and if cold gas is abundant what is the cause of the low rate 
of star formation? Focusing on these questions we investigate the 
effect of energetic particles on the cold gas <100 K in the filaments. 


(Johnstone 


3 MODEL 

The particle heating model used in this paper is essentially that de¬ 
scribed in F2009 for optically-thin emission from warm ionised gas 
adapted to include physics pertinent to very cold clouds. We do not 
consider here the ‘extra-heating’ case also considered in the afore¬ 
mentioned paper. The original model is motivated by spectra from 
extended emission line regions in BCGs which show spatially co¬ 
incident emission from both molecular and ionised gas, suggesting 
that the total emission is from clouds with a range of densities and 
temperatures. F2009 parametrise the contribution of the different 
phases to the total volume as a weighted sum of ‘cloudlets’ of dif¬ 
ferent densities but constant pressure. The weighting function, a, 
is taken to be a power-law in the cloud density. The pressure in the 
ionised filaments has been measured by [Heckman et al.[ ( |l989[ ) and 
is similar to the X-ray derived pressure of the ICM. 

[Ferland et al.[ ( [200^ use the Ha and H 2 12 fim emission line 
ratios to determine a for the Perseus cluster (a constant pressure 
of ~10®'^ cm“^ K is used for the Perseus cluster). They find an 
a — —0.35 can reproduce all emission lines from the warm ionised 
gas to within a factor of two (but see [Johnstone et al.[[2012[ for a 
discussion of the Ha flux in the Perseus cluster). This value of a 


implies the majority of the volume is filled by hot ionised gas with 
only a minor contribution from the more massive, cold clouds. Ha 
emission is predominantly from ~10^ K gas while the H 2 12 jim 
emission line is from gas ~10^ K. The authors point out that the 
data available at the time is not able to constrain a large reservoir 
of very low temperature clouds. However, extrapolation of a to low 
temperature clouds gives good agreement with mass estimates from 
CO emission lines though these may be optically thick ([Salome[ 
[etal.|2011[ |. Extrapolating this model of optically-thin lines to the 
very cold <100 K gas; the simple model predicts [O l]A63/im/ 
[C ll]A157/imratios of ~20. However, the observed line ratios are 
^1 ( Mittal et al.|2011 2012 Werner et al.|2013 2014). A myriad 
of physical processes could be responsible for this discrepancy in 
the observed ratio. 

The calculations in this paper were made using Cloudy C13 
( [Ferland et al.[[2013[ l and we provide our code as Appendix A of 
this paper. As in F2009, due to the longevity of the filaments, we 
neglect regions in which the cooling is thermally unstable. We cau¬ 
tion here that this may not be the case in the cold filaments and the 
line ratios in some lines are sensitive to this assumption (see Section 
[3.1[ ). [Chatzikos et ^ ( [2015[ ) studied the effects of non-equilibrium 
cooling but found it not to be important in the hot and warm gas. 

The models shown here are of a cloud where the farthest edge 
from us is not open, that is photons cannot escape out the far side of 
the cloud. While un-physical we choose to show these models due 
to simplicity of explaining the underlying physics of the emission 
lines. We have repeated all our calculations with an open-ended 
cloud, the only difference in results occurs in the last zone where 
the temperature of the open ended cloud drops as photons can es¬ 
cape. This drop in temperature has the effect of slightly decreas¬ 
ing the emissivity in the emission lines which remain optically-thin 
such as Ha. No significant effect is observed in the integrated line 
ratios. 

Cloudy contains a grain model which resolves the size dis¬ 
tribution of the grains and can calculate the grain properties, such 
as the charge distribution, temperature and grain opacities sepa¬ 
rately for grains of a given radius. The grain photoelectric heat¬ 
ing can then be determined for graphite, silicate and PAH features 
separately. A description of the Cloudy grain model is given in 
[van Hoof et al.[ ( |200T][2004[ l. We include dust grains with an ISM 
abunda nce and PAHs using the distribution from [Bakes & Tielens| 
( |l994| ) with an abundance of 10“^ ® carbon atoms per hydrogen 
atom in PAH features. The Jura rate is assumed for the formation 
of H 2 on grains and the dissociation of H 2 is calculated self consis¬ 
tently within Cloudy. 

In F2009 simulations are of a representative unit volume of 
gas. The actual filament is envisioned as being composed of many 
such parcels with a range of densities, temperatures, and ionizations 
and also with a range of column densities. For this approach to 
work the simulation of the unit volume is set to be optically thick 
in all resonance lines, as would happen if part of a much larger 
column density. CLOUDY calculates a full model of the hydrogen 
atom and its emission physics. We use this model in the predictions 
in this paper but assume that the Lyman lines are all optically thick. 
The effect of this optical depth is to stop fiuorescent pumping by 
the metagalactic background, and prevent Lyman lines produced 
in the recombination process from freely escaping. Rather, Lyman 
lines undergo a large number of scatterings and are converted into 
Baimer and other lines, are absorbed by dust, or escape from the 
cloud. 

The variable parameters in the particle heating model are the 
density of ionising particles measured relative to the Galactic cos- 
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Figure 2. The ionisation structure of hydrogen, oxygen and carbon. The solid contours indicate where the atomic species and the ionised/molecular species 
have the same abundance. The dashed contours where the gas is 90 per cent atomic and the dotted contours where the gas is 98 per cent atomic. 
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Figure 3. The log emissivity of the Ha (left), [O i] A63/4m (middle) and [C ii]A157/4m (right) emission lines for the simple model. The contour of constant 
gas pressure, relevant to the Perseus cluster (10®-^ cm“^ K) is drawn as a solid black line on each plot. The dashed black lines indicate the excitation 
temperature of the line from ground while the dotted black lines indicated the critical density of the emission line. In the optically-thin model we integrate the 
emissivity of each line for a given pressure over a weighted density distribution. 


mic ray background and the volume density of hydrogen in the gas, 
n. The energy density of the Galactic cosmic ray background is 
~1.8 eV cm“^ which corresponds to a pressure of ionising par¬ 
ticles of ~2xl0^ K cm“^. In the centre of the Perseus cluster 
[Sanders & Fabian| ( |2007| ) find an electron energy density ~1000 
times the Galactic background. We investigate the effect of ionis¬ 
ing energetic particles over a range of values corresponding to ~ 1 
to times the density of the Galactic cosmic ray background 
and gas densities of to hydrogen nuclei per cm cubed. 


3.1 Integration limits 

As mentioned in the previous section, the calculations in this paper 
neglect regions in which the cooling is thermally unstable. Thermal 
stability for a constant pressure gas is given by. 


'd{C -H)' 


'diC-HY 

Po 

'd{C-Hy 

dT 

p 

dT 

P ^0 

dp 


( 1 ) 

where C is the cooling rate, H is the heating rate, P is the pressure, 
p is the density and T is the temperature ( |Field|1965] see also Fig. 


13 of F2009). The left hand panel of Fig. shows the emissivity 
versus temperature profiles of Ha and many lines from cold gas 
for the zero turbulence case. The bottom panel show the cooling 
multiplied by the temperature. The thermally unstable regions of 
the grid are shown by the shaded regions and the line emissivity 
in these regions, or equivalently cloudlets with these temperatures, 
are not included in the integral. 


Additionally we use a lower and upper integration limit for 
the temperature of the cloudlets of the lowest temperature in the 
grid at that pressure and 10^ K respectively. The emission from 
cool and cold gas has very little emissivity beyond a few 10^ K 
so the choice of upper limit of the integral has negligible effect 
on the total emission in the lines. The lowest temperature, at each 
density, occurs where there is no additional ionising particle fiux. 
At this point the excitation in the gas is only by the background 
cosmic ray ionisation radi ation (we assume a m ean H° ionisation 
rate of 2x 10“^® from Indriolo et al.|2007| ), the CMB, and the 
ISM photoionising radiation field. We caution that the total emis¬ 
sion from some lines from very cold gas is sensitive to this lower 
limit cut-off. 
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Figure 4. Left: The temperature dependence of the emissivity of some FIR and sub-mm cooling lines. The regions where the gas is thermally unstable are 
shaded. These are determined by from the cooling xT curve in the bottom panel. Right: The line ratio, [O i]A63/rm / [C ii]A157/rm, calculated for an 
optically-thin cloud for different hydrogen densities and ionising particle densities. The ratio is very sensitive to variations in the cloud properties. The contour 
of constant total pressure, relevant to the Perseus cluster (10®-^ cm“^ K) is drawn as a solid black line. The ratio increases towards the top right-hand corner 
of the grid. This corresponds to an increase in density at constant temperature and can be understood as the critical densities of the [O i] lines are much larger 
than that of the [C ii] lines. 


3.2 Overview of physical and chemical properties of the grid 

Ionised and cold gas filaments in BCGs are situated in the hot 
ISM, surrounded by the X-ray emitting galaxy halo. The penetrat¬ 
ing nature of X-rays means they can input energy locally into cold 
and dense clouds creating partially ionised regions within predom¬ 
inantly atomic or molecular gas. The X-ray photons action is to; 
one, ionise; two, excite; or three, heat the gas. The contribution of 
each of these mechanisms depends on both the ionisation fraction 
(Xe = Tie In) of the gas and on its composition; whether the gas is 
mostly atomic or molecular. It is less sensitive to the shape of the 
ionising continuum, and the absolute values of the electron (rie) or 
hydrogen (n) densities. Additionally, electrons and ions with a high 
kinetic temperature exist in the hot halo which can also penetrate 
the filaments and ionise and excite the cold gas. 

In neutral gas most of the ionisation and excitation by the 
supra-thermal electrons will be of H° due to both its larger abun¬ 
dance and larger cross sections for interaction compared with other 
neutral elements or first ions. Ionisation and excitation of He° and 
He^ will be the next most important interactions. Primary and 
Auger ionisation will lead to secondary electrons which can col- 
lisionally ionise or excite other species in the gas. Excitation will 
generate line photons which can both interact further with the gas or 
may escape the cloud causing the gas to cool. Finally, direct heating 
of the gas is achieved through Coulomb collisions (see [Maloney,] 
[Hollenbach & Tielens|1996| for an in-depth description of X-ray ir¬ 
radiation of cold dense clouds). Here we are modelling the effect of 
collisional excitation with energetic particles which may arise from 
energetic particles in the hot halo interacting with the cold cloud or 
energetic electrons from X-ray irradiation of the cloud. 

The temperature, total pressure and ionisation fraction of the 
gas are shown in Fig.[2for the simple, optically-thin particle heat¬ 
ing model; no additional pressure from turbulent motions or mag¬ 


netic fields are included. In all three plots the thick black contour 
indicates a gas pressure of ~10®‘^ cm“^K, appropriate to the 
Perseus cluster of galaxies. As expected, as we increase the flux of 
ionising particles, for fixed density, the temperature and pressure 
increase and as does the ionisation fraction. For a fixed flux of ion¬ 
ising particles, increasing the density of the cloudlet corresponds 
to a decrease in the temperature and an increase in the pressure. As 
the temperature decreases so does the ionisation fraction. 

Where the ionisation fraction (Xe = Uejn) is greater than a 
few percent. Coulomb collisions with the thermal electron popu¬ 
lation dominate the interactions of the supra-thermal electrons and 
provide the majority of the heating in the cloud. When Xe falls 
below 1 per cent ionisations of neutral hydrogen dominate with ex¬ 
citations al so becoming increasingly important jShull & van Steen- 1 
|berg| 1985] [Maloney, Hollenbach & Tielens|1996| l. The ionisation 

structure of hydrogen, oxygen and carbon, for a range of ionising 
particle densities and hydrogen densities is shown in Fig. S H° per- 
sists over much of the parameter space with the gas becoming fully 
ionised at several tens of thousands of Kelvin and fully molecular 
at only a few tens. In atomic gas, when the energy of the electrons 
drops below the smallest excitation energy, all of the remaining en¬ 
ergy must go into heating the gas. Therefore, for a gas dominated 
by neutral hydrogen, once the supra-thermal electron energy falls 
below 10.2 eV, Coulomb heating is the only heating process. So a 
low ionisation fraction in atomic gas leads to a low heating effi¬ 
ciency ( [Shull & van Steenberg|1985[[Xu & McCray|199r| l. 

At high densities and low temperatures the gas transitions to a 
primarily molecular phase, though pockets of ionised and atomic 
regions will remain. Heating is more efficient in molecular gas 
and is dominated by collisional de-excitation of vibrationally ex¬ 
cited H 2 molecules and by H 2 dissociation by photons and ioni¬ 
sation by secondary electrons. Significant amounts of energy can 
also go into the excitation of the rotation-vibration bands of H 2 and 


(c) 0000 RAS, MNRAS 000, 000-000 





























6 R.E.A. Canning et al. 


into the dissociative electronic states ^Glassgold & Langer|1973| 
the temperature increases around a density 
of 10^'^ cm“^ . This is due to the increase in the H 2 gas fraction 
leading to a higher heating efficiency. 

Fig. shows the variation in the emissivity of the Ha, 
[O l]A63/imand [C ll] A 157/im emission lines respectively across 
the whole parameter space. The dashed line in each case indi¬ 
cates the excitation temperature of the line from the ground while 
the dotted black lines indicated the critical densities (for colli¬ 
sions with electrons) of the lines; ric? = 5 x 10^ cm“^ and 
ng "> = 3 X 10^ cm ^ . However, it should be noted electrons 
are not the only colliders in the gas. 

Collisional excitation is the dominant contributor to the line 
emission in all regions of the plot. However, above the excitation 
temperature excitation is by both thermal particles and secondary 
electrons whilst below the emissivity is dominated by excitation 
with secondary electrons. In our model Ha emission in the fila¬ 
ments arises from collisions with both populations (see upper panel 
of Fig. 15 of |Ferland et al.|206^ , and the geometry of the resulting 
clouds will be permeated by pockets of high and low temperature in 
a ‘Swiss cheese-like’ fashion. The degree of polarisation expected 
from the model is therefore complex (see also jSparks et al.|2014| |. 

The emissivity of the collisionally excited lines is essentially 
proportional to the UamfiT) where Ua and rih are the species col¬ 
liding and f{T) is the Boltzmann factor of the line for the case 
where collisions are occurring with thermal particles and is a con¬ 
stant for the case where the collisions are occurring with supra- 
thermal secondary electrons. Therefore, in the lower right hand cor¬ 
ner of the plots, where the gas is mostly molecular and the tempera¬ 
ture is too low for collisional excitation by thermal particles, at con¬ 
stant density the emissivity will increase by a more-or-less constant 
amount which is proportional to the number of secondary electrons 
available to collide with. The approximate fiattening of the emissiv¬ 
ity for constant ionising particle density with increasing hydrogen 
density is due to the high abundance of targets, such as H°, com¬ 
pared with the low abundance of ionising secondaries; increasing 
the density doesn’t affect the emissivity. The supra-thermal parti¬ 
cles prevent the gas from becoming fully molecular even in the 
densest regions. 

As the temperature of the gas increases and we move towards 
the upper left hand side of the plot collisions by thermal particles 
also become important and the line emissivity is less sensitive to 
the number of ionising particles. Significantly above the excita¬ 
tion temperature the emissivity remains essentially constant with 
increased gas temperature until the temperature is high enough that 
most atoms are ionised and the line emissivity drops. 

4 COMPLICATING THE MODEL 
4.1 Additional cold clouds 

The formalism presented in F2009 describes the cloud in terms of a 
weighted sum of very small ‘cloudlets’ of different density phases 
but at constant pressure. The weighting function is parametrised 
as a fraction of the total volume filled by each phase and an as¬ 
sumption is made that the function takes the form of a power-law 
in density. The total emissivity is found by integrating the density 
dependent emissivity over the cloud distribution. 

In reality we do not know the shape of the cloud distribution 
which may be highly complex; in practice a power law is chosen for 
simplicity. Physically, the variation of the power-law index (a) tells 


us the contribution to the total emissivity in a particular emission 
line, from gas of different temperatures, thus a relatively high a 
would show regions where dense, cold gas dominates over warmer 
less dense regions. As the cloud is assumed to maintain a constant 
gas pressure a simple relation exists for the temperature and total 
particle density. 

In the right hand panel of Fig.j^the ratio of [O l] A63/im over 
[C II] A157/im is predicted for a grid of ionising particle fluxes and 
cloud densities. Our assumption of constant pressure constrains us 
to lie along the black contour in the right hand panel of Fig. 
The hashed region shows where contours of constant pressure of 
10® cm“^ K to 10^ cm“® K lie. If a power law is not an appro¬ 
priate parametrisation for the coldest clouds then their predicted ra¬ 
tios will be different. For example, at a pressure of 10®'® cm“® K 
in gas with a temperature of T<25 K (n > 10®'^) [C ll]A157/im 
dominates over [O l] A63/im emission while at higher temperatures 
the contrary is true. 

Making the assumption that the power law is not an appropri¬ 
ate parametrisation for the coldest gaseous phases the most efficient 
way to produce a low ratio of [O l]A63/im over [C ll]A157/im 
would be by assuming additional cold clouds where the ratio is very 
low. The left hand panel of Fig. shows for a constant pressure 
cloud with 10®'® cm“^ K, the [O l]A63/im over [C ll]A157/im 
ratio is lowest (~ 0.05) for hydrogen densities of 10®-® cm-'^ 
and where the ionising particle density is of order 10^'® times the 
Galactic cosmic ray background. 

Assuming a cloud distribution with the power law function 
described in |Ferland et al.| ( [T009| ) and an additional delta function 
where the ratio of [O l]A63/im over [C ll]A157/im is lowest we 
require the emissivity of this additional cold cloud to be a factor 
of 1.8 X 10^ times that expected for a 1 cm^ cloud with the same 
properties, in order for the [O l]A63/im over [C ll]A157/im ratio to 
be ~ 1. However, the addition of more cold clouds has a significant 
effect on the other line ratios. Importantly Ha and H 2 12.28 /xm are 
significantly stronger for high density, low temperature clouds than 
[C ll]A157/im by factors of ~ 100 and ~20,000 times respectively, 
which is not observed in the filaments (see Fig.[^. 


4.2 Optically thick lines 

The majority of the lines of interest in F2009 are optically-thin thus 
validating the use of a unit volume of gas in the predictions of the 
line ratios. However, as previously mentioned, some lines in the 
wavelength range of Herschel and ALMA may become optically- 
thick at reasonable column densities and as such require special 
attention. 

The intensity of radiation from a cloud decreases along its path 
due to scattering and absorption and increases due to spontaneous 
and stimulated emission. For an optically-thin plasma, all radiation 
generated within the plasma is able to leave it, spontaneous emis¬ 
sion is all that is required and the matter is not in thermal equilib¬ 
rium with the radiation. Hence an optically-thin source will radiate 
below the black body limit. In an optically-thick system photons 
are absorbed and re-emitted many times before leaving the cloud 
and the cloud is considered to be at a quasi-equilibrium temper¬ 
ature where emission from the cloud in all directions will be the 
same. This sets an upper limit to the intensity obtained at a specific 
frequency from a thermal source at a specific temperature and is 
the black body limit. The luminosity in the line, per steradian, at 
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Figure 5. The ratio of the emissivity of Ha to that of [C ii]A157/im (left) and the ratio of the emissivity of H 2 12.28 fim to that of [C ii]A157/im (right) 
for optically thin clouds with the same conditions as in Fig.|4| Observationally Ha/[C ii]A157/im ~1.3 (Mittal et al. 1 20 12|| Werner et al.|2014) and H 2 12fi 
m/[C ii]A157Aim- 0.025 (F2009). 


the black body limit can be estimated by, 

Lune = TbfvAVTrR^ =(’^] - fvAVnR^ ( 2 ) 

\kBj -1 

where Af/ is the velocity width of the line in km , due to mi¬ 
cro turbulence, fv is the velocity filling factor, and Tb is the line’s 
brightness temperature in K. For [O l]A63/im and [C ll]A157/im 
we cannot make the simplifying assumption that hu kT. 

[Mittal et al.| ( |2011| ) have observed [C ll]A157/im and [O 
l]A63/im line emission using Herschel and measure line widths 
of a few hundred kms“^. The flux expected at the black body 
limit is therefore W m~^ for [C ll] and 

W for [O l] (with a filling factor of 1). The measured flux of 
[C ll]A157/im is comparable to that of the black body limit while 
the measured flux of [O l]A63/im is a factor of one hundred lower. 

This is a very conservative limit as the real velocity width of 
the line is likely to be a lot smaller. [Salome et ^ < |2008b) measure 
line widths of the CO J = 2 — 1 line in NGC 1275 in the Perseus 
cluster with the 2.5” beam of the Plateau de Bure Interferometer, to 
be as low as ~30 kms“^. The smallest Herschel angular resolu¬ 
tion is 5”, while the cool ionised gas filaments are determined, with 
the Hubble Space Telescope, to have an upper limit of only 0.2 arc- 
sec (70 pc) wide in Perseus. It is therefore very likely, within this 
large aperture, that the velocity width is dominated by the macro 
turbulence of many intertwined filaments, not the micro turbulence 
of the emission line. 

To explore the effect of high column on the predictions of 
the particle heating model we will first make the simplifying as¬ 
sumption that the cloud is at one initial temperature and model 
the emergent line intensities of important IR cooling lines for 
clouds of different column density. In this manner we are study¬ 
ing conditions across a finite column density cloud with surface 
properties corresponding to a single point in Fig. We show 
in Fig. I^the emission line ratios of [O l]A63/im/[C ll]A157/im 
and [O l]A63/im/[0 l]A145/im for four such models with surface 
temperatures of 15 K, 33 K, 60 K, and 100 K at a pressure of 
~ 10®'^ cm“^ K, the emissivity profiles of many cold gas lines 
are shown in Fig. jBlj 


Fig. shows that for a constant pressure cloud, even at low 
initial temperature, [O l]A63/im/[C ll]A157/im > 1 once the lines 
are optically-thick. This is expected from the black body limit, as¬ 
suming the emission lines come from the same gas, given in Equ. 
[^ If both lines are saturated we expect [O l] A63/im/[C ll] A157/im 
~3, therefore under equilibrium conditions and assuming no ad¬ 
ditional pressure support in the gas, an optically-thick gas is not 
sufficient to explain the observed line ratios of [O l]A63/im and 
[C ll]A157/im in the filaments. 

The [O l]A63/imand [O l] A 145/im emission lines should be 
produced from the same clouds so this ratio has less degenera¬ 
cies than the [O l]A63/im/[C ll]A157/im ratio and can be a more 
sensitive test of the optical depth in the lines. In the optically 
thick limit, for a reasonable range of temperatures. Fig. [^ shows 
the ratio of [O l]A63/im/[C ll]A157/im should be lower than has 
currently been observed (observations are indicated by the grey 
region). We note that in many cases [O l]A145/im is too weak 
to be detected and as such the upper limit to the grey region is 
likely to increase with deeper observations. In addition, in the 
no-turbulence case for lines with temperature <100 K the [O 
I]63/im line becomes optically-thick around a column density of 
~ 10^^ cm“^, corresponding to a high extinction of Ay =5.6 
(assuming Nh=1.8x lO^^Ay from Predehl & Schmitt||l995| l The 
[C ll]A157/imand [O l]A145/imlines become optically thick at 
columns which are prohibitively large (see also jLiseau, Justtanontj 
|& Tielens|2006j ). 


4.3 T\irbulence 


In the previous sections, and in the F2009 model, no turbulence 
was included in the lines. We currently only have upper limits on 
the turbulent motions of a single cloud in the extended filaments 
of BCGs as the beam size is large compared with the filament 
threads. Optical and NIR observations of the filaments typically 
measure line width s, of gas at >^1,00 0 and 1 0,000 K, to be between 
50 — 200kms“^ ( Lim et aLl|2012j ) while Salome et al. ( |2008bt 
have measured line widths of ~30 km , in the cold gas, using 
the 2.5” beam of the Plateau de Bure Interferometer in the extended 
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Figure 7. The ratios of key lines overlaid with contours of gas pressure in the cases where none, 2 km and 10 km turbulence is present in the gas. 
Left: [O I] A63)um/ [C ii] A157)um. Middle: Ha / [C ii] A157/im. Right: H 2 12.28 fim / [C ii] A157/im. The measured fluxes are the integral of the cloudlets, 
indicated by the pressure contour. Observationally the fluxes are [O i]A63yum/ [C ii]A157/im< 1, Ha/[C ii]A157/im ~1.3 jMittal et al.| 2012[ [Werner et al.| 
and H 2 I2ii m/[C ii]A157/im- 0.025 (F2009). 


filaments of NGC 1275. The 2.5 arcsec beam is ~12 times the up¬ 
per limits on filament thread width determined by the Hubble space 
telescope ( [Fabian et al.|2008| l. All these observations currently suf¬ 
fer from confusion of filament threads. 

In Galactic giant molecular clouds line widths are typically 
measured as a couple of km in the molecular gas with a mea¬ 
surements of a few km s~^ in the cold and warm neutral mediums 
( jHennebelle & Chabrier||2013j ). Under 100 K the sound speed is 
already less than 2 km while in < 30 K gas the sound speed 
is below 1 kms“^ so molecular clouds may be turbulent in the 
molecu lar and possibly cold neutral phases. [Kritsuk & Norman] 
( |201l| ) show that scaling relations found between the ^^CO(J=1-0) 
line-width and size and mass and size of molecular clouds ( |Lar-| 
|son|198l] l can be interpreted as a signature of supersonic motion. 
Other supporting evidence for supersonic motions come from the 
measurement of a log normal distribution of the column density 
( |Vazquez-Semadeni||1994| l. However, while turbulence predicts a 
log normal distribution it is not necessary that a log normal distri¬ 
bution implies turbulence ( jTassis et al.|2010| ). 

Micro-turbulence increases the line width which suppresses 
the optical depths and alters the importance of shielding and pump¬ 
ing of lines. Self-shielding becomes less important while the ability 
to absorb a larger part of the continuum increases the importance 
of fiorescence. Our initial grids, of the emissivities of the emission 
lines, shown in Fig.j^are not affected by the addition of turbulence 
as in this simple model all lines are assumed optically-thin. We note 
here our cloudy models do not include shock heating of the gas. 

As well as altering the chemical balance, the introduction of 
microturbulence affects the pressure balance of the cloud. Cloudy 
includes the turbulence as a velocity in kms“^, Uturb- The addi¬ 
tional energy density is therefore Pturb — ^pUturb^ where p is the 
gas density and F is a constant which accounts for how ordered 
the turbulent velocity field is. We assume F = 3, appropriate for 
isotropic turbulent motions but note that the filaments are likely 
threaded with magnetic fields which would infiuence the turbulent 
motions of ions ( jHeiles & Crutcher|2005| l. 

The effect of the additional pressure term significantly alters 
the predicted line ratios, from cold gas, for reasonable values of 
Uturb- The gas pressure contour required to keep the filament in 
pressure equilibrium with the surrounding hot gas is reduced as 


shown in Fig. Essentially, when the turbulent pressure is large 
the gas cools at constant density. The turbulent pressure is density 
dependent and as such the largest deviation from the no-turbulence 
curve, for a constant turbulent velocity at all temperatures, is in the 
densest gas phases. As can be seen in Fig. lathis will alter the ratios 
of [O I] A63/im/ [C ll] A157/im, pushing them to lower ratios while 
leaving those of higher temperature lines unchanged. 

Turbulent dissipation would also lead to heating of the gas. 
We can estimate the contribution of the turbulent heating to the gas 
line luminosity using, 

Pturb — 2 ^ cvgs s , (3) 

where Mtot is the gas mass, vturb is the turbulent velocity and I 
is the injection scale of the turbulence. We assume the very con¬ 
servative values of Mtot ~ 10®M© and vturb ~ 10kms“^. 
The scale is the largest unknown, not least because individual fil¬ 
ament threads have not been resolved, and so we assume here 
/ ~ 1/4 HST filament width limit ~ 17 pc. These values lead to 
Lturb ~ 5.5 X 10^^ ergs s“^. In the horseshoe filament chosen as 
a typically example filament by F2009 in NGC 1275, the authors 
find a Ha luminosity of 7 x 10®® ergs s“^. The [C ll]A157/im/ Ha 
ratio is ^0.8 so L^cii] ^ 5.6 x 10®® ergs s“^. The turbulent heat¬ 
ing is therefore only one percent of the line luminosity. A turbulent 
velocity of ~45 kms“^ is required to reach the [C ll]A157/im 
line luminosity, or alternatively I ~ 0.17 pc. However, this still 
neglects the luminosity of H 2 , CO and other important lines. 


4.4 Magnetic fields 


The presence of strong magnetic fields in the filaments of BCGs has 
been inferred from optical observations of the geometry and widths 
of single threads of extended filaments ( [Fabian et al.|2008t , from 
radio observations ( Taylor et al.[2006[ ) and implied from density ar¬ 
guments ( [Werner et al.[2013| . However, magnetic fields have never 
been directly detected in the filaments. If these strong fields exist 
they will increase the energy density, Pmag — also decreasing 
the contribution to the total pressure budget from gas pressure. 

In our implementation B is a constant and not dependent on 
density. This may not be the case but a full treatment is beyond 
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Figure 8. The ratios of key lines overlaid with contours of gas pressure in the cases where none, 10 juG and 100 juG fields are present in the gas. Left: 
[O i]A63)um/ [C ii]X157/j,m. Middle; Ha / [C ii]X157/im. Right; H 2 12.28 /im / [C ii]X157/im. The measured fluxes are the integral of the cloudlets, 
indicated by the pressure contour. Observationally the fluxes are [O i]A63yum/ [C ii]A157/im< 1, Ha/[C ii]A157/im ~1.3 jMittal et al.| 2012[ [Werner et al.| 
and H 2 12/1 m/[C ii]A157/im- 0.025 (F2009). 


the scope of the paper. As the B field is not density dependent the 
additional pressure is a constant in all cloudlets and as such the 
required gas pressure is less at all densities. This is shown graph¬ 
ically by the contours on Fig. Importantly, as in the case of ad¬ 
ditional turbulent support, the required density of the coldest gas is 
less dense than in the fiducial model which decreases the predicted 
ratios of [O l]A63/im over [C ll]A157/im. However, this imple¬ 
mentation will effect all line ratios; the optical and IR lines from 
the 10,000 K—1,000 K gas as well as the FIR and sub-mm lines 
from much colder gas. 

4.5 Words of caution 

The models presented above explore the effect of the relaxation 
of certain assumptions, namely, the assumption of a power-law in 
density, optically-thin gas and no turbulent or magnetic pressure 
support to the gas. However, we relax these assumptions one-by- 
one. In nature we may expect more than one of these effects to be 
present in the gas which will clearly introduce degeneracies but in 
the interests of presenting simple testable models we do not explore 
the entire parameter space. It is therefore important to consider this 
when devising observational strategies to test the model. We dis¬ 
cuss the models and what observations may be robust tests of these 
models in the next section. 


5 PREDICTIONS AND DISCUSSION 

Extended cool/cold gaseous filaments are observed in BCGs yet, 
in most cases, apparently have very low star formation efficiencies. 
The question of what may be preventing the gas from cooling - 
what is its excitation/heating mechanism - is an important one as it 
goes to the heart of how much can cooling hot X-ray gas affect late 
stage galaxy evolution. In this paper we hope to provide testable 
predictions for the ‘particle’ excitation mechanism and to elucidate 
the importance of the surrounding hot gas in inhibiting the cold¬ 
est/densest gas from forming and hence preventing star formation. 

Making the assumption that the model put forward in F2009 
and explored further here is the dominant excitation mechanism 
of the ionised and cold gas phases we can make predictions for 


further high spatial resolution cold gas observations of these fila¬ 
ments. The high spatial resolution of ALMA and HST is impor¬ 
tant as deriving gas masses and measuring the turbulent veloc¬ 
ity of the filaments is highly dependent on the assumed gas frac¬ 
tion within the beam and on whether we are observing one or 
many filament threads. In Section we look at the effects, on the 
[O l]A63/im and [C ll]A157/imline ratios in the model, of relaxing 
one-by-one some of the assumptions of F2009. Many mechanisms 
may be playing a part in setting the line ratios. However, applying 
Occam’s razor we wish to single out a possible dominant mecha¬ 
nism. 

Observations have found the ratios of [C ll]A157/im/Ha is 
approximately 0.8, [O l]A63/im/[C ll]A157/im <1 and where it 
has been possible to measure the typically weak [O l]A145/im 
emission line the ratio of [O l]A63/im/[0 l]A145/im~5-30 
|tal et al.||2011| |2012[ [Werner et al.||2013| |2014| ). As discussed in 
Sections [4~T] the addition of a sink of very cold clouds at tempera¬ 
tures ~10 K cannot explain the FIR line ratios while maintaining 
the ratios observed in the optical and NIR if the gas is all at constant 
pressure with the surrounding X-ray gas. [C ll]A157/im becomes 
optically thick at a greater column than [O l]A63/im (see Apendix 
[B] ) and the high [O l]A63/im/[0 l]A145/im line ratios also rule out 
these emission lines being optically thick and suggest they come 
from the warm 50—100 K gas (see Fig.|^. Masing and foreground 
absorption in the [O l] lines can change these ratios but would tend 
to lower rather than increase the ratio (see discussion in |Liseau,| 
jJusttanont & Tielens|2006| l. 

Assuming the clouds maintain constant pressure with the hot 
ICM the most promising scenario would be to include an addi¬ 
tional, density dependent, pressure term. Increasing the pressure 
support in this density dependent fashion, such as through an 
increase in the turbulence or density dependent magnetic fields 
pushes the ratio of both [O l]A63/im/Ha and [C ll]A157/im/Ha 
to lower ratios for a range of pressures and powerlaw indexes. An 
additional pressure support such as this would leave the pressure 
contribution from the hot gas essentially the same while requiring 
less thermal pressure from the dense gas, changing the line ratios 
that are predicted from cold, dense clouds while not altering the line 
ratios from the 10,000 K, ‘warm’, gas. If turbulence contributes all 
the additional pressure support then 2-10 km turbulence is re- 
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Figure 9. Left: If the gas is optically thin then at reasonable temperatures and densities the [O i]A63)Um / [C ii]A157/im ratio traces fairly well the gas 
pressure. Therefore at constant temperature (or density) we can use the ratio as a density (or temperature) diagnostic. Here we plot the emission lines ratio 
against the density for clouds of different temperatures. The grey band indicates the detected line ratios in some BCGs and giant elliptical galaxies. Right: A 
[O i]A63)Um / [O i] A145)um versus [O i]A63)Um / [C ii]A157/im diagnostic diagram for a cloud at 60 K. 


quired to reach the observed line ratios of approximately unity. This 
additional support decreases the density of the cold gas and could 
therefore also explain the longevity of the filaments. 

Other indications exist that thermal pressure support alone 
is unlikely to be able to keep the filaments in BCGs stable for 
long periods of time. The periodicity observed in the star forma¬ 
tion regions of the filaments of NGC 1275 would require turbu¬ 
lence of ~17 kms“^ if all the additional pressure support came 
from turbulence ( [Canning et al.|2014] ). However, we note that the 
linear structure and long lifetimes suggest some pressure support 
likely comes from magnetic fields. Additionally, in radio-galaxies. 


quiescent, out-flowing cold clouds are often observed (e.g. 

Nes^ 

jvadba et al.||2010| [Combes et al.[[2013| [Alatalo et al.[[2014 

1 , al- 


beit at greater velocities than filaments from typical ‘mechanical 
AGN feedback’ in BCGs. jGuillard et aLl ( |2014| l show that turbulent 
energy seeded by the AGN feedback could be responsible for the 
bright [C ll]A157/im emission lines observed in these sources. 

Whilst the far infrared lines observed by Herschel offer impor¬ 
tant diagnostics of the gas properties, estimating the total mass in 
the very cold gas still requires an extrapolation of the model to the 
coldest emission lines < 20 K. Most current observations of CO 
have observed the low J levels which are likely to be optically thick 
and therefore unable to constrain the gas masses. Future observa¬ 
tions targeting optically thin CO lines will help answer the question 
of how much mass is concealed in the cold gas in massive galaxies 
and therefore the efficiencies of their star formation rates. 


5.1 Model predictions 

The model we are advocating in this paper is one where the exci¬ 
tation is dominated by collisions with energetic particles creating 
pockets of partially ionised regions interspersed with molecular re¬ 
gions in the filaments in a ‘swiss-cheese-like’ fashion. The gas is 


supported also by magnetic fields or modest turbulence, which it¬ 
self could be seeded by the impinging particles, which add pressure 
support to the gas, lowing the density and preventing immediate 
collapse into stars. In this section we outline tests of this model and 
diagnostics of the conditions of the gas. 


5.1.1 [O l]X63iim/[0 l]X14:5iim and 
[O l7A63/im//C ll7A157/im 

Our models show, for optically thin emission, in low temper¬ 
ature gas <150 K the [O l]A63/im/[C ll]A157/im ratio traces 
relatively well the pressure. So for fixed temperature the ra¬ 
tio can be used as a density diagnostic. In the left hand panel 
of Fig. 1^ we provide a plot of the optically thin line ratios 
of [O l]A63/im/[C ll]A157/im against density for gas with tem¬ 
peratures ranging between 15 K and 100 K. The current val¬ 
ues of the ratio from Herschel observations are indicated by the 
grey region. Taking the gas to be ~60 K, consistent with the 
[O l]A63/im/[0 l]A145/im in the optically thin limit, in the right 
hand panel we plot the [O l]A63/im/[0 l]A145/im against the 
[O l]A63/im/[C ll]A157/im ratios and indicate the densities. Cur¬ 
rent observations, where at least two of these lines are detected, 
are overplotted. For gas under these conditions the current observa¬ 
tions indicate the density should be a few 10^ cm“^. It is important 
to note these ratios are for the total galaxy fiuxes which includes 
the denser and in some cases star forming interiors as well as the 
extended filaments. 


5.1.2 High J CO emission lines and high density gas tracers 

Observations of the high J (J=3-2 or higher) CO rotation ladder 
with ALMA will shed further light on the excitation mechanism 
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Figure 6. The predicted [O i]X63fj,m / [C ii]A157yL4m (top) and 
[Oi]A63/im/ [O i]A145/im (bottom) line ratios versus column density 
for cloudlets at a constant pressure of ~ 10®'^ cm“^ K, but with 
varying surface temperatures. At high optical depths the [O i]A63)Um / 
[C ii]A157yum ratio tends to a few, this is larger than is observed. Obser¬ 
vational [O i]A63)Um / [C ii]A157/im ratios in the filaments are typically 
less < 1. In the bottom plot the grey region shows the currently detected 
values for the [O i]A63)Um/ [O i]A145)Um line ratio, we note that in many 
cases the [O i]A145/im emission line is too weak to detect. These line ra¬ 
tios favour optically thin gas at ‘warm’ temperatures (50-100 K). 


and the importance of both turbulence and magnetic fields. Observ¬ 
ing high J transitions are important for three key reasons. First, the 
high J transitions are insensitive to the low temperature integration 
limits (see Fig.[^; second, they are less likely to be optically-thick 
than the low J CO lines; and third, the shape of the high J rota¬ 
tion ladder will enable a sensitive determination of the gas tem¬ 
perature; we might expect the gas temperature to differ to the dust 
temperature if the gas is excited by energetic particles as is seen in 
the molecular clouds in the nuclear disk of the Galaxy (e.g.|Yusef-| 
IZadeh, Wardle & Roy|2007l l. 

The ratios of these high J CO lines to Ha is also a sensi¬ 
tive test of the presence of gas with densities greater than 10^ 
(see Fig. which allows limits to be placed on the levels of ad¬ 
ditional pressure support allowed as the assumption of total pres¬ 
sure equilibrium therefore enforces constant density cooling at very 
low temperatures. High density lines like those of HCN and HCO^ 


(n > 10^ cm“^) should not be observed in the extended non-star 
forming filaments if the additional pressure sup port is equivalent to 
10 kms“^ turbulence. ISalome et al.|(|2008a|) have detected HCN 


(3-2) in the central regions of NGC 1275 but so far no detection of 
very dense gas in the extended regions of filaments has been ob¬ 
served. Fig. |C2| shows the expected emissivities of CO J(3-2), CO 
J(5-4), HCN and HCO^ for a range of ionising particle fluxes and 
gas densities, the emissivity in these lines is highly density depen¬ 
dent. Constraining the density of the gas allows us to identify the 
level of additional pressure support and whether the clouds remain 
in pressure equilibrium in the coldest < 20 K gas. Additionally, 
high spatial resolution spectroscopy with ALMA will enable us to 
measure the turbulent support in the cold gas. Putting these together 
we can also elucidate the level of magnetic support required to the 
total pressure budget. 


5.1.3 Low J emission lines 


Optical depth estimates from the ^^CO J(l-O) emission lines sug- 
gest this line is optically thick in the fi laments of BCGs ( |Salome| 
& Combes||2003 Salome et al.||2011 note that ^^CO and C^^O 
are likely to be optically thin in the filaments and as such may be 
good diagnostics of the model). Ratios of these lines with other 
detected emission lines are thus not easily compared to predictions 
from our model. However, these lines are likely thermalised and 
can be used as sensitive thermometers for the gas as the excitation 
temperature will be equal to the kinetic temperature. Combined 
with ratios of high J emission lines of CO which are sensitive 
tracers of the density (see Fig. \C2) constraints can be put on the 
physical conditions of the coldest gas. 


We have argued above that collisional excitation with the hot 
particles in the surrounding X-ray gas may inhibit the dense con¬ 
ditions required for star formation in the extended filaments. How¬ 
ever, we know of a few cases in which some filaments have dis¬ 
rupted into stars (e.g . [McDonald & Veilleux|r2009[ |0’Dea et al.| 
|2010| [Canning et al.||2014| |. In regions in which young stars are 

observed obviously emission lines dominated by excitation from 
young stellar sources will be observed. We might also expect the 
turbulence to be higher in these regions due to driving from young 
stellar winds hence the diagnostics suggested above are not rele¬ 
vant to regions which are actively forming stars and should not be 
applied to star forming knots in filaments. 

However, some information may still be gleamed from study¬ 
ing these regions. Additional pressure support of the gas, either 
by turbulence or by magnetic fields should alter the Jeans mass. 
Collisional excitation with energetic particles from the surround¬ 
ing gas heats the gas and increases its ionisation fraction which 
will increase its coupling to magnetic fields, which should slow the 
timescale for gravitational collapse. We speculate that this excita¬ 
tion mechanism could therefore lead to a top heavy IMF in the fil¬ 
aments explaining the very massive star clusters observed in some 
outer filaments ( [Canning et al.|2014] l. 


6 CONCLUSIONS 

We have extended the particle heating mechanism first presented 
in [Ferland et al.[ ( [2009[ ) to examine the effects of abundant cold 
clouds, optically-thick gas, turbulence and magnetic fields on the 
bright line ratios from cold, <100 K, gas. We explore these models 
in the context of extended cool and cold gaseous nebulae observed 
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in some massive galaxies many of which are devoid of star forma¬ 
tion. We wish to be clear that we are not considering the centrally 
peaked emission and associated star forming regions observed in 
many massive galaxies where it is clear many process are occur¬ 
ring. We focus instead on ‘clean’ non-star forming, extended nebu¬ 
lae. 

We suggest the simplest explanation for the discrepancy be¬ 
tween the predicted [C ll] and [O l] line ratios may be that there is a 
small amount of additional pressure support in the cold gas from ei¬ 
ther turbulence (2-10 km ) and/or density dependent magnetic 
fields and we present predictions for line ratios and diagnostics of 
the gas temperatures and densities. We suggest that turbulence may 
be driven by the infiux of energetic particles impinging on the cold 
gas, creating a higher level of ionisation in the gas and decreasing 
the gas density which could also inhibit star formation in these fil¬ 
aments. The majority of filaments would remain long lived unless, 
the gas pressure drops, star formation is either triggered by exter¬ 
nal perturbations, such as disturbances from rising AGN bubbles or 
galaxy interactions or the heat from such external perturbations is 
strong enough to eradicate the filaments. 
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APPENDIX A: CLOUDY CODE 

c 

c 

table HM05 z=0 
extinguish by 21, leakage = 0 
cmb redshift 0 
c 
c 

atom H2 levels large 

atom H-like Lyman pumping off 

abundances he =-1.022 li =-10.268 be =-20.000 b =-10.051 c =-3.523 n 
=-4.155 

continue o =-3.398 f =-20.000 ne =-4.222 na =-6.523 mg =-5.523 al =-6.699 
continue si =-5.398 p =-6.796 s =-5.000 cl =-7.000 ar =-5.523 k =-7.959 
continue ca =-7.699 sc =-20.000 ti =-9.237 v = -10.000 cr = -8.000 mn 
=-7.638 

continue fe = -5.523 co =-20.000 ni =-7.000 cu =-8.824 zn =-7.6990 no 

grains 

grains ism 

grains pah 

set pah constant -4.6 

set H2 Jura rate 

case B 

c 

c 

cosmic rays background -1 vary 
grid from 0 to 7 in 0.1 dex steps 
hden 4 vary 
grid 0 6 0.1 
c 
c 

stop zone 1 
set dr 0 

stop temperature off 
turbulence = 2 km/s 


APPENDIX B: OPTICAL DEPTHS 


In Rig. |B1| we shown the emission line emissivity divided by the 
cloud column density against the cloud column density for two 
models. The first has a hydrogen density of 10^ ^ cm“^ and a par¬ 
ticle density of 10^ times that of the Galactic background cosmic 
ray energy density (the Galactic background is here taken as 1.8 eV 
cm“^) which results in a cloud surface temperature of ~15 K. We 
choose this temperature and density model as an optically-thin gas 
under these conditions has [O l]A63/im over [C ll]A157/im less 
than one. The second has a hydrogen density of 10^ ° cm“^ and 
a particle density of 10^'^ times that of the Galactic background 
cosmic ray energy density leading to a higher temperature, ~33 
K. This is close to the temperature at which the emission expected 
from both lines is equal for an optically-thin gas (see right hand 
panel of Rig.|^. Both models have a pressure of ~ 10®'^ cm“^ K. 

Ror an optically-thin line the emissivity will increases in pro¬ 
portion to the column, this produces horizontal lines on the plots. 
The column density at which the column begins to increase faster 
than the emissivity indicates the column at which the line has be¬ 
come optically-thick and line photons are collisionally dexcited 
following multiple scatterings. The [O I]63/im line in the shown 
model (upper left plot), becomes optically-thick around a column 
density of ~ 10^^ cm“^, correspondi ng to a high extinction o f 


Av =5.6, assuming Nh=L8x10^^Av ( Predehl & Schmitt] 1995 1. 
The [O l]A145/imand [C ll] A 157/im emission lines become opti¬ 
cally thick at columns greater than ~ 10^^'^ cm“^ which are very 
rare and so these lines are likely optically thin. 
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In the right hand panels we show the emissivities in the CO 
rotation ladder as a function of column density. The low J lines of 
CO have large optical depths in gas with low extinction Ay < 1 and 
are most likely optically thick. These lines may not be constraining 
to test the various excitation models of the filaments but can be 
sensitive thermometers of the gas temperature. 


APPENDIX C: EMISSIVITIES 

For completeness in Fig. |C1| and |C2| we show the optically thin 
emissivities of other key strong lines from our model and the ratios 
for the neutral carbon and oxygen lines. In the cold gas the optically 
thin [O I] line ratios are a relatively faithful indicator of the gas 
temperature while the [C l] ratios are more sensitive to the particle 
ionising flux. The lines shown in Fig. |C2| are good indicators of the 
gas density. 
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Figure Bl. The predicted line emissivity divided by the column density plotted against the column density of the clouds for clouds at constant pressure. The 
surface temperature of the cloud in the upper plot is T~33 K while the lower plot has a lower surface temperature of T~15 K. 
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Figure Cl. Predicted log emmisivities of the [O i]A145iim, [C i]A369/im, and [C i]A609/im transitions and the neutral carbon and neutral oxygen line 
ratios for an optically thin gas. The countour indicates a gas pressure of 10®-^ cm“^ K . 
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Figure C2. The predicted log emissivities of the CO J(3-2) and J(5-4) transitions and those of lines of HCN and HCO+ overlaid with contours of the gas 
pressure as shown in Fig. [7] All these lines offer sensitive diagnostics of the gas density and the CO ladder transitions can also offer sensitive temperature 
measurements. If a high degree of additional pressure support is present in the filaments than the densest gas should not be. Observations of the strength of 
these density sensitive lines in regions in which no star formation is occurring can constrain the additional pressure support in the gas. 
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